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Abstract 


Pulsed  jets  have  been  shown  to  penetrate  into  a  cross  flow  more  than  a  corresponding 
steady  state  jet.  The  increased  penetration  depends  to  some  extent  on  the  formation  of 
distinct  ring  vortices  which  in  turn  transport  themselves  into  the  cross-flow  by  the  self 
induced  velocity.  This  work  has  been  concentrated  in  understanding  the  physics  of  the 
increased  penetration.  Experiments  have  indicated  the  formation  of  vortex  rings  as  a  jet 
was  pulsed  even  when  the  propane  fuel  was  combusting.  The  dynamics  and  trajectories 
of  vortex  rings  formed  by  pulsation  of  a  jet  in  a  uniform  cross-flow  were  studied  and 
detailed  measurements  were  made  using  flow  visualization  techniques  including  laser  in¬ 
duced  fluorescence  and  hot-film  anemometry.  Measurements  indicated  that  vortex  rings 
were  fully-formed  at  a  distance  of  three  times  the  jet  exit  diameter.  A  numerical  exper¬ 
iment  based  on  a  Lagrangian,  grid-free,  three-dimensional  vortex  method  was  performed 
and  general  agreement  with  the  experiments  was  concluded.  At  low  frequencies,  inter¬ 
action  between  adjacent  vortex  rings  was  negligible  because  of  relatively  large  separation 
distances  and  under  these  conditions,  each  ring  behaves  like  a  single  discrete  vortex  ring.  It 
was  observed  that  single  vortex  rings  moved  into  the  flow  sometimes  tilted  up  to  about  30°, 
depending  on  the  ring’s  strength.  Tilting  into  the  flow  further  contributed  to  the  increased 
penetration.  Numerical  simulation  showed  that  the  tilting  is  due  to  the  combined  effect 
of  viscosity  and  the  cross-flow.  Single  pulses  into  a  supersonic  cross-flow  also  indicated 
the  formation  of  vortex  rings.  Difficulties  were  experienced  in  accomplishing  continuous 
pulsed  injection  into  a  supersonic  cross-flow. 
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1,  Introduction 


Efficient  mixing  of  two  fluids  frequently  plays  a  major  role  in  fluid  engineering  prob¬ 
lem  affecting  combustion,  propulsion,  petrochemical,  food  engineering,  and  environmental 
applications.  In  these  applications,  a  fluid  is  frequently  introduced  in  the  form  of  a  jet 
into  a  main  flow.  The  spread  of  a  fluid  jet  into  an  ambient  stream  occurs  at  a  fixed  rate 
which  is  to  slow  under  certain  conditions.  As  a  result,  mixing  enhancement  techniques  are 
frequently  introduced  into  transverse  jet  flows  to  increase  the  mixing  rate  between  the  jet 
and  the  surrounding  fluid.  Eric  and  Roshko  (1989)  investigated  different  types  of  vortical 
structures  which  exist  in  the  near  field  of  the  transverse  jet.  They  showed  that  the  wake 
vorticity  comes  from  the  cross-flow  boundary  layer,  and  not  from  the  jet  fluid.  The  result  of 
their  experiment  implies  that  the  wake  may  not  significantly  contribute  to  mixing.  Broad- 
well  and  Breidenthal  (1984)  noted  that  the  transverse  jet  appears  to  be  a  better  mixer  of 
jet  and  cross-flow  than  a  free  jet.  An  important  technique  for  enhancing  the  mixing  of 
axial  jet  flows  has  been  to  introduce  different  periodic  disturbances  into  the  mixing  process 
to  increase  the  rate  of  mixing.  Pulsation  of  axial  jets  for  mixing  enhancement  was  first 
introduced  through  the  development  of  the  pulse  jet  engine  (Lockwood,  1963).  He  noted 
the  generation  of  vortex  rings  by  pulsating  flow.  Vandsburger  at  al.  (1986)  also  observed 
the  formation  of  large-scale  vortex  rings  in  the  flame  zone  by  exciting  a  laminar  flame  at 
a  frequency  of  16  Hz.  This  problem  has  received  increased  attention  in  terms  of  mixing 
improvements  (Sarohia  and  Bernal,  1981;  Reynolds  and  Bouchard,  1981;  and  Bremhorst 
and  Harch,  1979). 

As  another  efficient  means  to  enhance  the  mixing  process  through  deeper  penetration 
and  increased  entrainment  rate,  Wu  et  al.  (1988)  introduced  a  transverse  injection  with 
forced  pulsation  into  a  uniform  cross-flow  with  unsteady  effects.  Introducing  unsteadiness 
in  the  jet  of  such  flows  generally  resulted  in  flow  patterns  and  features  which  were  totally 
different  from  those  of  a  steady  jet  and  which  could  be  of  considerable  value  in  improving 
the  mixing  rate  and  increasing  the  entrainment  and  diffusion  of  jets  in  cross-flow.  The 
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objective  of  the  present  study  was  to  investigate  the  increased  mixing  due  to  the  pulsation 
of  the  jet  in  cross-flow.  At  low  frequencies,  vortex  rings  are  formed  penetrating  into  the 
cross-flow  with  relatively  high  speeds.  In  order  to  accomplish  the  above  purpose,  we  need 
to  understand  the  basic  mechanism  for  vortex  ring  formation,  dynamics  and  penetration 
into  the  cross-flow. 

Vortex  ring  formation  and  motion  have  been  the  focus  of  research  by  many  investi¬ 
gators  for  over  a  century  ever  since  Helmholtz  (Thomson,  1883)  formulated  the  general 
vortex  theorems.  The  motion  of  vortex  rings  is  one  of  the  interesting  fluid  phenomena 
with  concentrated  vorticity.  The  fact  that  a  vortex  ring  can  transport  fluid  and  momen¬ 
tum  over  a  certain  distance  has  been  one  of  the  motivations  for  vortex  ring  research.  The 
early  works  on  the  vortex  ring  phenomena  are  well  described  in  the  paper  by  Guhler  and 
Sallet  (1979).  These  early  studies  (Thomson,  1883  and  Hicks,  1885)  are  based  on  the  fol¬ 
lowing  simplifying  assumptions;  the  ratio  of  core  radius  to  ring  radius,  a/jR,  is  very  small 
and  the  vorticity  is  concentrated  in  the  vortex  core.  Magarvey  and  Maclatchy  (1964)  used 
photography  to  extensively  investigate  the  formation,  structure  and  disintegration  of  vor¬ 
tex  rings  by  flow  visualization  with  smoke  in  air.  The  propagation  property  of  vortex  rings 
in  air  has  been  investigated  by  using  stroboscopic  photography  and  hot-wire  anemometry 
by  Johnson  (1970).  The  study  was  devoted  to  accurately  measuring  the  decay  of  real 
vortex  rings  in  quiescent  air  and  comparing  calculated  ring  behavior  with  the  various  the¬ 
oretical  predictions.  Sallet  and  Widmayer  (1974)  measured  the  characteristic  features  of 
laminar  and  turbulent  vortex  rings  with  hot-wire  anemometers.  These  measurements  were 
comprised  of  the  translational  velocity,  ring  diameter,  core  diameter,  and  circulation  at 
several  points  downstream  from  the  orifice  exit.  The  thorough  discussion  of  the  structure 
and  stability  of  viscous  vortex  rings  was  qualitatively  performed  in  a  series  of  papers  by 
Maxworthy  (1972,  1974,  1977).  Using  flov/  visualization  and  laser  doppler  velocimetry  he 
investigated  the  formation  process,  vorticity,  the  flow  fields,  the  translational  velocity,  and 
stability  of  vortex  rings  over  a  wide  range  of  Reynolds  numbers.  The  temporal  development 
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of  the  rolling-up  process  of  vortex  rings  was  investigated  in  detail  using  flow  visualization 
and  laser  doppler  velociraeter  by  Bidden  (1979).  He  examined  the  variation  of  velocity, 
diameter,  and  circulation  of  vortex  rings  which  move  between  the  point  of  formation  and 
a  location  three  jet  diameters  downstream. 

Many  studies  have  been  performed  based  upon  vortex  methods  to  simulate  a  vorticity- 
dominated  flow  of  an  incompressible,  two-dimensional  flow  since  Rosenhead  (1931).  Even 
though  the  method  has  some  inherent  difficulty  that  has  limited  the  application  of  the 
method,  it  is  generally  known  that  the  rotationed  region  of  fluid  motion  is  well  simulated 
by  a  number  of  discrete  vortex  filaments  in  an  in  viscid  flow.  A  vortex  ring  is  a  typical 
example  of  vorticity- dominated  flow  with  a  small  compact  region  of  vorticity.  There  have 
been  a  number  of  simulations  of  two-dimensional  axisymmetric  vortex  rings  performed;  but, 
some  studies  were  done  on  three-dimensional  vortex  ring  configuration  with  variations  of 
size,  position,  and  shape  considered  together.  A  series  of  vortex  rings  can  be  simulated  by 
the  three-dimensional  discrete  vortex  method  and  the  time  dependent  motion  of  each  of  the 
discretized  vortex  elements  can  be  traced  throughout  the  flow  field  with  this  computational 
methodology.  Detailed  analysis  and  review  of  the  discrete  vortex  method  have  been  given 
by  Leonard  (1980,  1985),  Ashurst  and  Meiburg  (1988)  and  Sarpkaya  (1989)  in  modeling 
jets,  wakes,  boundary  layers,  vortex  rings,  and  other  vorticity- dominated  flow. 

To  the  authors’  knowledge,  there  has  not  been  a  systematic  study  of  the  dynamics  of 
vortex  rings  in  cross-flows.  Therefore,  this  study  was  divided  into  two  parts.  First,  through 
an  experimental  investigation,  the  vortex  rings  produced  by  pulsed  transverse  injection 
into  a  uniform  cross-flow  were  studied  using  flow  visualization  techniques  and  hot-film 
anemometry.  Second,  a  numerical  simulation  for  three-dimensional  vortex  rings  in  cross- 
flows  was  carried  out  using  the  Lagrangian,  grid-free,  vortex  method  to  help  understand 
and  to  compliment  the  experimental  results. 
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2,  Experimental  Procedures 


2.1  Description  of  the  Water  Tunnel  Facility 

The  water  tunnel,  located  at  The  University  of  Tennessee  Space  Institute,  is  a  closed- 
return  type  and  the  dimensions  of  the  test  section  are  30.5  cm  high,  45.7  cm  wide  and  150 
cm  long  (12  in  x  18  in  x  60  in).  The  tunnel  speed  can  be  varied  from  1.5  cm/sec  to  60 
cm/sec.  The  test  section  walls  are  constructed  of  1.27  cm  thick  Plexiglas  for  versatility 
in  observing  and  photographing  the  flow  fleld.  The  measured  turbulence  intensity  in  the 
test  section  was  measured  as  0.1%.  The  nozzle  between  the  stagnation  chamber  and  test 
section  is  a  bell  mouth  type,  which  gradually  and  continuously  changes  from  a  circular  to 
rectangular  cross-section  with  a  contraction  ratio  of  13.5.  The  tunnel  velocity  field  was 
calibrated  by  using  a  cylindrical  hot-film  probe  to  ensure  uniform  cross-flow.  The  probe 
was  calibrated  by  using  a  specially  constructed  water  channel  (Chang,  1991). 

2.2  Experimental  Setup  for  Continuous  Pulsation 

Figure  1  shows  the  schematic  of  the  experimental  setup  in  the  water  tunnel.  The  jet 
was  installed  on  a  false  floor  of  the  test  section  which  provided  a  known  boundary  layer 
at  the  jet  location.  The  false  floor  was  a  surface-finished  plate  with  an  elliptic  leading 
edge  whose  aspect  ratio  is  2:1  and  it  was  located  19  mm  above  the  tunnel  floor  and  was 
1500  mm  long.  The  jet  exit  (dj  =  12.7  mm)  was  located  190  mm  behind  the  leading  edge. 
This  allowed  the  laminar  boundary  layer  thickness  at  the  jet  to  be  on  the  order  of  one  jet 
diameter.  The  Reynolds  number  based  on  the  jet  diameter  and  steady  jet  velocity  was 
Rtdj  =  1.87  X  10^.  The  tunnel  floor  and  one  side  wall  were  marked,  using  a  square  grid 
pattern,  with  lines  spaced  at  equal  distance  to  be  used  with  the  flow  visualization  data 
collection  of  vortex  rings  size  and  trajectories. 
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Pulsation  of  the  full  jet  flow  was  accomplished  by  incorporating  a  solenoid  valve  inline 
in  the  water  supply  tube.  This  valve  was  driven  by  a  frequency  generator  so  that  any  fre¬ 
quency  pulsation  could  be  investigated.  The  valve  did  not  respond  uniformly  at  frequencies 
of  20  Hz  and  above.  The  solenoid  valve  was  approximately  20  diameters  upstream  of  the 
jet  exit;  thus,  ensuring  a  uniform  pipe  flow  at  the  jet  exit  (Chang,  1991). 

The  motion  of  vortex  rings  was  studied  for  several  different  area-averaged  jet  exit 
velocities  to  cross-flow  velocity  ratios,  U j/Uoo  —  1-5  to  6.7.  The  jet-equivalent  velocity 
was  14.88  cm/sec  which  is  a  velocity  without  pulsation  in  the  pipe.  The  pulsed  jet  was 
maintained  at  a  flow  rate  of  18.85  cm^/sec  (0.3  GPM),  therefore,  the  variations  in  UjIUoo 
were  obtained  by  the  changes  in  the  cross-flow  velocity  Uoo-  The  jet-equivalent  velocity  in 
the  pipe  should  be  distinguished  from  the  jet  exit  velocity,  (U e),  caused  by  the  pulsation 
in  the  same  pipe.  This  jet  exit  velocity  was  close  to  the  initial  translational  velocity  of  the 
vortex  ring. 

Typical  velocity  traces,  at  the  jet  exit,  at  different  frequencies  are  presented  in  Figure 
2.  The  original  pulsed  velocity  profile  generated  by  a  function  generator  was  a  square 
wave.  The  low  frequencies  (e.g.,  1  and  2  Hz)  pulsations  make  the  flow  more  like  a  periodic 
train  of  impulsively  started  rings  in  a  cross-flow  rather  than  an  excited  jet  in  a  cross-flow. 
The  velocity  trace  characteristics  have  a  significant  effect  on  actual  vortex  ring  formation 
as  does  the  piston  velocity  profile  in  a  single  pulsation  by  the  piston-cylinder  assembly. 
As  the  frequency  increases  the  maximum  velocity  (impulse)  is  considerably  reduced.  This 
means  that  at  high  frequencies  the  circulation  of  the  ring  is  considerably  decreased  while 
a  steady  flow  rate  remained  constant. 

2.3  Hot- Film  Anemometry  and  Data  Acquisition 

Hot-film  anemometry  was  used  to  measure  the  velocity  field  and  determine  charac¬ 
teristic  parameters  of  vortex  rings  in  cross-flow.  A  single  hot-film  probe  was  utilized  with 
anemometers  operated  in  the  constant-temperature  mode  to  measure  mean  velocities  in 
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the  three-dimensional  flow  field.  The  Constant  Temperature  Anemometry  (CTA)  system 
is  comprised  of  a  normal  TSI  hot-film  probe,  a  TSI  model  IFA  100  Intelligent  Flow  An¬ 
alyzer,  a  Scientific  Solution’s  Lab  Master  16  channel  A/D  converter,  a  1  kHz  low  band 
pass  filter  and  a  IBM  PC  with  20  megabyte  hard  disk.  An  oscilloscope  was  connected  to 
the  control  box  to  observe  the  voltage  traces  obtained  from  measurements  within  the  flow 
field. 

Figure  3  represents  a  schematic  diagram  of  the  data  acquisition  system.  The  voltage 
signals  from  the  hot-film  probe  were  sent  into  the  Intelligent  Flow  Analyzer,  (IFA),  which 
was  connected  to  the  Lab  Master  A/D  converter.  The  signals  passed  through  the  1  kHz 
band  pass  filter  to  eliminate  noise.  The  signals  were  digitized  and  were  sent  to  a  VAX 
11/785  computer  for  storage  and  processing  to  obtain  the  velocity  components  from  the 
calibration  results. 

Probe  calibration  is  described  in  Chang  (1991)  with  a  discussion  of  directional  charac¬ 
teristics  of  a  hot-film  probe.  Also,  general  expression  for  the  effective  cooling  velocity  of  a 
probe  in  the  three-dimensional  flow  field  is  derived  for  the  interpretation  of  the  anemometer 
output  signal. 

2.^  Flow  Visualization  Techniques 

Two  types  of  flow  visualization  techniques  were  utilized  for  water  tunnel  experiments 
to  complement  each  other  for  a  complete  description  and  analysis  of  the  flow  details.  The 
techniques  unveiled  the  global  flow  features  of  vortex  rings  in  cross-flow.  The  flow  was 
visualized  using  a  common  dye  technique.  The  dye  was  a  mixture  of  commercial  food 
coloring,  alcohol,  and  evaporated  milk  with  a  specific  gravity  of  one.  Internal  dye  ports 
were  built  into  the  jet  exit  region  to  allow  flow  visualization  with  minimal  disturbance 
into  the  local  flow  and,  therefore,  the  vortex  rings  were  colored.  Both  a  camcorder  and 
conventional  35  mm  camera  with  1/1000  sec  frame  setting  were  used  to  take  pictures. 
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A  laser  light  sheet  technique  was  applied  to  provide  a  sectional  view  of  the  vortex 
rings.  An  Argon  ion  laser  light  source  with  a  maximum  power  output  of  2  watts  was 
used  for  the  laser  sheet  visualization.  A  45°  inclined  mirror  and  a  4  mm  cylindrical  lens 
were  employed  to  reflect  and  produce  a  plane  laser  light  sheet.  Rhodamine  6G  fluorescent 
dye  was  supplied  in  small  concentrations  through  four  small  holes  around  the  nozzle  to 
coordinate  the  pulsing  fluid.  Sectional  views  of  the  vortex  ring  were  obtained  to  analyze 
the  detailed  structure  of  the  ring  and  its  wake. 

2.5  Measurement  Uncertainties 

The  experimental  uncertainties  associated  with  measurements  of  the  mean  jet  velocity 
(U j),  sensor  operating  temperature  (T,),  corrected  effective  cooling  velocity  (U eff,c)^  time- 
averaged  mean  velocity  components  (U  and  W^),  the  measurement  of  circulation  (r^) 
and  ring  size  (R)  were  estimated  for  a  typical  vortex  ring  at  a  frequency  of  1  Hz  the 
Z  —  3dj.  Following  are  the  experimental  uncertainty  estimates  of  the  above  properties 
actually  measured.  AUj  :  ±  11.0%  (this  is  a  maximum  possible  error,  based  on  usual 
error  analysis,  and  in  reality  could  be  much  less),  ATg  :  i:  0.15%,  AU eff,c  '  i  0.57%,  AU 
:  ±  2.82%,  AW  :  ±  0.60%,  Ar„i  :  ±  12.42%,  AR  :  ±  4.57%.  These  estimates  are  valid 
for  all  vortex  ring  cases  with  the  same  odds.  The  uncertainties  in  the  calculated  value  was 
evaluated  using  the  standard  uncertainty  equation  (Kline  and  McClintock,  1953)  based  on 
the  uncertainties  in  the  primary  measurements  (Chang,  1991). 
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3.  Experimental  Results 


3.1  Motion  of  Vortex  Rings  in  Crossflow 

Figure  4  shows  selected  still  photographs  which  reveal  several  characteristic  vortex 
ring  features  observed  in  the  present  study  such  as  the  size,  tilting,  and  trajectory  of  the 
vortex  rings  while  moving  in  the  cross-flow.  From  these  observation  it  was  evident  that 
vortex  rings,  for  the  range  of  parameters  in  this  study,  move  as  blobs  of  fluid  with  little 
change  in  their  shape.  The  tilted  shapes  and  wake  structures  of  the  vortex  rings  were 
identified  using  laser-induced  fluorescence  (LIF)  flow  visualization.  The  core  of  the  vortex 
ring  is  defined  as  the  region  where  the  velocity  varies  linearly  with  the  distance  from  its 
center.  A  wake  region,  where  a  small  amount  of  the  ring  fluid,  therefore  vorticity,  is  left, 
also  exists  behind  the  ring. 

As  the  pulsation  frequency  increased  in  this  experiment,  two  separate  physical  phe¬ 
nomena  were  seen  to  influence  the  vortex  rings  and  their  trajectories.  First,  the  circulation 
strength  within  each  ring  was  reduced,  which  directly  resulted  in  the  reduction  of  its  in¬ 
duced  translational  velocity.  This  enabled  the  cross-flow  momentum  to  overcome  each  ring 
easily  and  deflect  it  into  the  downstream.  Second,  the  distance  between  the  sequential  rings 
reduced  as  the  frequency  increased.  This  resulted  in  interactions  between  the  rings,  such 
that  the  vortex  rings  quickly  coalesced. 

In  this  case,  when  two  sequential  vortex  rings  are  close  to  each  other,  the  vortex  ring 
diameter  of  the  forward  vortex  ring  increases  with  decreasing  vortex  ring  velocity;  the 
diameter  of  the  rearward  vortex  ring  decreases  while  its  vortex  ring  velocity  increases. 
Since  this  interaction  is  present  during  the  initial  development  of  the  individual  rings,  the 
rings  are  not  fully-formed  as  they  are  turned  into  the  downstream  direction.  The  jet  flow 
development  into  the  cross-flow  including  the  jet  trajectory,  at  high  pulsing  frequency  of 
larger  than  15  Hz  is  approaching  that  of  the  steady  jet  case. 


Figure  5  shows  the  trajectories  of  the  vortex  rings  at  various  frequencies  for  different 
jets  to  cross-flow  velocity  ratios.  These  vortex  ring  trajectories  exhibit  two  fundamental 
results.  First,  the  penetration  level  at  low  frequencies  is  much  higher  than  at  high  frequen¬ 
cies  for  the  constant  jet  to  cross-flow  velocity  ratios;  second,  for  a  given  pulsing  frequency 
the  penetration  height  decreases  as  the  jet  to  the  cross-flow  velocity  ratio  decreases. 

3.2  Tilting  Mechanism  for  a  Vortex  Ring  in  Crossflow 

Vortex  rings  are  generated  by  the  periodic  opening  and  closing  of  a  valve  in  the  supply 
line  into  a  cross-flow  environment.  Vortex  rings  are  formed  at  the  exit  of  the  line,  if  the 
opening  and  closing  takes  place  in  short  time  pulses.  At  low  frequencies,  rings  are  distinctly 
separated  from  each  other  with  minimum  interactions.  Vortex  rings  produced  in  this 
manner  move  as  unit  blobs  for  several  diameters  into  the  flow.  The  motion  is  interesting 
in  that  it  involves  tilting  of  the  ring  into  the  cross-flow.  The  tilting  process  was  initially 
believed  to  be  related  to  the  interferences  which  might  exist  between  sequential  rings. 
However,  experiments  showed  that  even  for  a  single  vortex  ring  there  was  an  appreciable 
amount  of  tilting  present  as  it  moved  from  the  exit  into  the  cross-flow.  Figure  4(f)  shows 
the  motion  of  a  single  vortex  ring  into  a  uniform  cross-flow.  The  tilting  is  most  pronounced 
for  rings  with  small  circulation,  generated  with  a  short  stroke.  At  higher  frequencies  the 
rings  are  weakened  and  are  spaced  very  close  to  each  other  and  the  tilt  is  increased,  thus, 
the  interaction  between  adjacent  rings  could  also  affect  the  tilting.  Pulsation  frequency 
affects  circulation,  size,  and  translational  velocity  of  the  rings,  therefore,  for  each  pulsation 
frequency  the  vortex  ring  will  be  formed  with  a  different  range  of  flow  variables. 

The  tilting  mechanism  of  a  vortex  ring  in  cross-flow  may  be  explained  using  a  qual¬ 
itative  approach.  At  low  frequencies  they  are  assumed  to  be  non  interfering:  therefore, 
each  vortex  ring  may  be  analyzed  individually.  Ting  and  Tung  (1965)  used  a  spinning 
solid  cylinder  as  a  model  for  the  dynamics  of  a  fluid  vortex  core  with  uniform  vorticity. 
They  used  the  Kutta-Joukowski  theorem  to  calculate  equivalent  forces  on  the  vortex  core. 
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However,  Caruthers  (1991)  has  suggested  that  a  Rankine  vortex  model  is  kinematically 
inconsistent  with  the  Kutta-Joukowski  result  and  that  a  singular  layer  of  vorticity  on  the 
surface  of  the  circular  cylinder  should  be  added  to  the  model  in  this  case.  Superimposing 
a  constant  vorticity  within  the  cylinder  (Rankine  vortex)  serves  to  transport  the  surface 
layer  vorticity  as  a  solid  cylinder  rotation  of  the  initial  distribution  and  its  potential  field. 
He  suggested  that  the  implied  lift  on  the  fluid  within  the  circular  cylinder  obtained,  based 
on  this  assumption,  was  only  one  half  of  the  Kutta-Joukowski  result  for  a  solid  body.  Even 
though  there  still  exists  a  controversy  about  the  magnitude  of  the  forces  and  the  tilting 
mechanism,  the  acting  force  is  proportional  to  the  Kutta-Joukowski  result  with  the  real¬ 
ization  that  the  actual  forces  depend  on  the  details  of  the  actual  core  and  layer  vorticity 
distribution.  A  non  uniform  interaction  with  the  cross-flow  leads  to  the  tilting  of  the  ring 
as  it  is  formed  near  the  jet  exit.  The  tilting  then  results  in  the  upstream  travel  of  the  ring 
as  it  moves  into  the  cross-flow,  as  shown  in  Figure  4(f).  This  figure  represents  a  single 
ring  moving  into  the  cross-flow  of  the  tunnel.  There  are  no  other  vortex  rings  in  the  flow 
and  the  tunnel  flow  is  uniform  except  for  a  relatively  thin  boundary  layer  on  the  tunnel 
floor.  The  tilting  and  subsequent  upstream  motion  of  the  vortex  ring  is  evidenced  from 
the  trailing  dye. 

3.3  Quantitative  Measurements 

A  large  amount  of  the  research  into  vortex  rings  is  based  upon  flow  visualization,  and 
are  therefore,  to  some  extent,  descriptive  and  subjective.  There  have  been  relatively  few 
quantitative  studies  to  measure  the  detailed  structure  of  the  ring.  The  circulation,  vorticity 
distribution,  core  size,  and  translational  velocity  can  be  determined  throughout  the  ring 
from  measurements  of  the  velocity  components.  In  order  to  better  understand  some  of 
the  above-mentioned  specific  features,  it  was  necessary  to  investigate  the  quantitative 
characteristics.  In  this  section,  the  results  of  hot-film  measurements  are  discussed. 
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To  illustrate  the  present  quantitative  measurements,  a  typical  vortex  ring  in  the  cross- 
flow  was  selected  with  a  velocity  ratio,  UjfUoo  =  3.0,  and  a  frequency,  /  =  1  Hz  where 
no  interaction  between  neighboring  vortex  rings  exists.  The  velocity  fleld  {W  and  of 
the  ring  in  cross-flow  was  measured  at  three  different  locations  (Z  =  3Dj,  7Dj  and  lOdj). 
The  mean  pulsating  velocity  (u)  at  /  =  1  Hz  was  12.63  cm/sec  which  corresponds  to  the 
velocity  ratio,  u/Uoo  =  2.4  in  this  case.  A  single  hot-film  probe  was  used  to  measure  two 
velocity  components,  W  and  [/,  at  a  fixed  point  as  the  vortex  ring  passes  by  the  measuring 
point.  Two  velocity  components  are  measured  by  rotating  the  plane  of  the  prong  through 
90°.  The  first  probe  position  selected  was  parallel  to  the  cross-flow  (9  =  0°)  and  the  other 
probe  position  was  perpendicular  to  the  cross-flow  (9  =  90°)  in  the  Y-Z  plane. 

The  tilting  angle  of  the  vortex  ring  was  determined  by  measuring  the  time  difference 
between  the  peak  velocities  at  the  front  inner  edge  and  rear  inner  edge  of  the  ring.  The 
tilted  angle  measured  by  this  method  was  about  3°  at  Z  =  Zdj.  The  center  position  of  the 
ring  was  X  =  OAOdj  at  Z  =  3dj,  which  means  the  ring  moved  sightly  toward  the  stream 
wise  direction  due  to  the  cross-flow.  Some  typical  velocity-time  traces  along  the  radial 
positions  at  Z  =  3dj  are  presented  in  Figure  6.  The  velocity-time  traces  were  obtained  by 
using  a  ^  =  0°  hot-film  probe.  Figure  6(a)  shows  a  typical  velocity-time  trace  when  the 
center  of  the  vortex  ring  passes  by  the  hot-film  probe.  The  dip  is  caused  by  the  rotational 
motion  of  the  ring  core.  Figures  6(b)  and  (c)  indicate  typical  velocity  traces  at  the  inner 
edge  of  the  ring  core  and  the  center  of  the  vortex  ring,  respectively. 

Tests  were  repeated  from  four  to  eight  times  at  each  of  three  downstream  locations  to 
ensure  the  repeatability  of  the  results  and  to  obtain  the  averaged  data  in  the  analysis.  Hot- 
film  measurements  at  Z  =  lOdj  were  much  more  difficult  to  do  than  at  Z  =  3dj  due  to  the 
various  ring  trajectory  from  run  to  run.  It  has  already  been  noted  that  some  data  scatter 
was  observed  relatively  far  downstream  in  the  flow  visualization  results.  It  appeared  that 
the  data  scatter  might  be  due  mainly  to  a  slight  non-uniformity  in  pulsations  generated  by 
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the  combination  of  the  solenoid  valve  and  the  frequency  generator.  Therefore,  great  care 
was  especially  taken  at  these  point  measurements  to  minimize  their  effect. 

Figure  7  represents  some  typical  non-dimensional  velocity  profiles  W /Ut  and  U /Ut 
versus  the  non-dimensional  coordinate  Z/D  (D  is  a  ring  diameter)  which  corresponds  to 
the  time  that  a  vortex  ring  passes  through  the  probe  at  Z  =  Zdj.  As  mentioned  above, 
the  ring  was  surveyed  along  nine  radial  positions.  Here,  the  distribution  of  each  velocity 
component  is  transformed  to  the  steady  coordinate  system  by  converting  t-coordinate  to 
Z-coordinate  using  Z  =  Ut  ■’t-,  where  Ut  is  the  instantaneous  translational  velocity  of  the 
ring  at  a  measuring  point.  Data  was  collected  at  a  rate  of  600  samples  per  second.  The 
translational  velocity  Ut  and  the  ring  diameter  D  at  the  measuring  point  were  selected 
as  characteristic  parameters.  The  assumption  of  quasi-steady  motion  implies  that  the 
properties  of  a  vortex  ring  are  not  significantly  changed  by  the  cross-flow  and/or  wakes 
during  the  time  that  a  vortex  ring  passes  by  a  measuring  point.  Since  the  hot-film  probe 
cannot  detect  the  negative  velocity  component,  the  negative  velocity  profiles  in  Figure  7 
are  resolved  by  flow  visualization. 

In  order  to  determine  the  translational  velocity,  core  diameter,  ring  diameter,  and 
circulation,  several  different  methods  can  be  adopted.  Locations  of  vortex  ring  core  centers 
were  very  carefully  determined  from  traces  like  those  in  Figure  6(a).  Near  the  core  center, 
the  probe  was  moved  in  small  increments  of  0.0156  cm  to  be  able  to  detect  the  core 
center  position  as  accurately  as  possible.  For  example,  the  core  diameter  at  Z  =  Sdj  was 
repeatable  to  within  a  maximum  ±4%  (Chang,  1991).  Then,  the  inner  edge  and  outer 
edge  of  the  core  were  determined  from  the  inner  peak  velocity  point.  The  core  radius 
between  the  position  of  peak  velocity  and  core  center  was  decided  by  assuming  that  the 
core  is  exactly  circular  in  shape  by  having  circular  streamlines.  It  was  confirmed  that 
the  core  diameter  obtained  from  the  above  method  is  equivalent  to  the  distance  between 
peak  velocities.  The  peak  velocity  at  the  inner  edge  Uu,  is  the  sum  of  rotational  velocity 
Urot  and  the  translational  velocity  Ut-  Also,  the  outer  edge  velocity,  Uoei  is  the  difference 
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between  Urot  and  Ut-  Therefore,  the  translational  velocity  and  rotational  velocity  of  the 
core  can  be  calculated  as  follows: 

UT=\{Uie  +  Uoe)  (1) 

Urot  =  liUie  -  Uoe)  (2) 

The  effective  core  circulation  can  be  calculated  from  the  core  rotational  velocity. 

Fc  =  t^rot  •  StTO  (3) 

The  total  measured  circulation  was  evaluated  by  integrating  the  velocity  profiles  W  as 
a  function  of  Z  along  the  vortex  ring  centerline.  Hence,  we  can  calculate  the  amount 
of  vorticity  contained  within  the  core.  In  this  case,  Tc/Tm  was  estimated  to  be  about 
65.5%  at  downstream  locations,  where  Fm  is  the  total  measured  circulation  obtained  from 
velocity  profiles  W. 

Table  1  presents  a  summary  of  quantitative  data  obtained  from  hot-film  experiments 
at  three  different  downstream  locations  for  /  =  1  Hz  pulsation.  The  Reynolds  number 
based  upon  the  ideal  circulation,  Rer,  ,  was  3970  in  this  case.  F,  indicates  the  initial  ideal 
circulation  calculated  using  the  simplified  slug  model  which  is  a  predicted  value.  Fj  was 
computed  by  using  the  following  relation  with  the  exit  velocity  profile  of  Figure  2. 

Fi  =  ^^  u‘^{t)di  (4) 

where  u{t)  is  the  mean  pulsating  velocity  of  the  fluid  ejected  through  the  exit.  The  amount 
of  measured  circulation  almost  linearly  decreased  as  the  ring  moved  downstream. 

We  need  to  discuss  an  important  point  concerning  the  relationship  between  measured 
core  circulation  and  measured  total  circulation.  Also,  it  is  significant  to  highlight  the 
discrepancy  between  the  measured  circulation  and  ideal  circulation  using  a  simple  slug 
model.  Maxworthy  (1977)  noted  that  only  about  50%  of  the  total  circulation  is  contained 
within  the  core  at  high  Rcm  based  on  the  exit  velocity  and  exit  diameter.  Contrasts 
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between  the  present  study  and  the  theoretical  result  of  Sallet  and  Widmayer  (1974)  who 
obtained  65.5%  and  71.5%  of  concentration  within  the  core,  respectively.  Maxworthy 
(1977)  also  found  that  the  total  of  a  measured  amount  of  vorticity  that  appears  downstream 
is  very  close  to  that  predicted  by  a  simple  slug  model  at  high  Reynolds  number  (between 
10^  and  10^).  However,  Bidden  (1977),  using  Laser  Doppler  Velocimeter  data,  observed 
that  Bi  underestimates  the  total  measured  circulation  by  30%  to  60%  at  low  Reynolds 
numbers  until  Rcm  =  6.6  x  10^.  In  this  case,  the  piston  velocity  profiles  were  utilized  to 
calculate  B,  .  This  also  indicates  that  the  simple  slug  model  cannot  be  a  good  parameter  to 
describe  the  formation  process.  Finally,  it  can  be  seen  that  the  above  relationship  between 
core  and  measured  circulations  depends  on  a  given  piston  motion  and/or  pulsed  frequency 
and  exit  geometry.  In  Didden’s  experimental  work,  the  measured  circulation  is  much  larger 
than  the  ideal  circulation  using  the  simple  slug  model.  Bidden  (1979)  explained  in  detail 
the  factors  that  can  increase  the  circulation  above  the  ideal  value  as  follows;  one  is  due  to 
a  much  larger  velocity  than  the  exit  velocity,  near  the  sharp  corner  at  the  nozzle  exit  for 
short  times;  the  other  is  due  to  the  accelerated  core  flow  for  constant  piston  velocity  by  a 
thicker  internal  boundary  layer.  In  the  present  study,  the  total  measured  circulation  was 
very  close  to  the  ideal  circulation  using  the  slug  model.  At  Z  =  10  dj,  it  deviated  only 
15.7%  from  the  ideal  quantity. 

The  tilting  angle  slightly  increased  as  the  ring  moved  downstream,  however,  it  was 
observed  that  the  tilting  angle  of  a  ring  tended  to  be  restored  to  the  horizontal  orientation 
farther  downstream  (near  the  free  surface)  even  though  it  was  not  quantitatively  confirmed 
by  a  hot-film  probe  experiment.  The  average  tilting  angle,  for  the  range  of  flow  parameters 
in  this  study,  was  estimated  to  be  approximately  3°  to  4°  during  the  vortex  ring  movement 
in  the  cross-flow.  Higher  tilting  angles  of  up  to  30°  were  observed  for  slower  flow  rates. 
The  expression  for  translational  velocity  by  in  viscid  flow  theory  is  expressed  as  follows: 
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which  is  comparable  to  the  quantity  from  vortex  ring  measurements,  a  non-dimensional 
velocity  was  introduced  as 

(6) 

i  M 

The  above  equation  is  equivalent  to  ln{8Rla)  -  1/4  in  the  in  viscid  equation.  Under  the 
present  experimental  range,  the  deviations  from  the  real  value  was  reasonably  small.  This 
makes  sure  that  the  results  using  equation  (5)  for  >1  =  1/4  for  a  viscous  vortex  ring  is 
accurate  within  30%. 

3.4  Measurement  of  a  Fully-Formed  Vortex  Ring 

In  order  to  quantitatively  investigate  the  distance  from  the  exit  a  fully-formed  vortex 
ring  is  produced,  traces  of  point-wise  velocity  versus  time  histories  were  obtained  at  twelve 
different  locations  (from  the  exit  to  Z  =  4dj)  along  the  jet  exit  centerline.  The  data  shown 
in  Figure  8  indicate  that  the  ring  is  basically  fully-formed  at  Z  =  8Dj  downstream  from 
the  exit.  Up  to  Z  =  2dj,  the  velocity  traces  look  like  jet  flow  movement  without  forming 
a  real  vortex  ring.  The  velocity  travels  between  2dj  and  Zdj  are  considered  in  this  study 
as  showing  the  formation  processes  of  a  vortex  ring.  When  a  fully-formed  vortex  ring  is 
passed  through  the  hot-film  probe,  the  velocity  profile  became  symmetric  about  a  peak 
velocity,  e.g.,  after  Z  =  3  dj,  as  shown  in  Figure  8.  After  the  ring  passed  downstream  of 
Z  =  3  dj  ,  no  change  in  the  configuration  of  its  continuing  formation  process  was  observed. 
Whitehead  (1968)  numerically  investigated  the  vortex  ring  formation  process  up  to  five  exit 
diameters.  He  suggested  that  2dj  is  enough  distance  to  obtain  a  fully-formed  vortex  ring. 
Sallet  and  Widmayer  (1974)  observed  by  photographic  study  that  the  rings  are  formed 
within  about  2.5dj  from  the  vortex  generator.  They  also  noted  that  the  early  development 
of  fully-formed  vortex  rings  seem  to  be  independent  not  only  of  the  piston  speed  but  also 
holds  true  for  vortices  formed  with  different  piston  velocity  profiles. 

Table  2  presents  some  comparative  data  on  vortex  ring  formation.  In  Table  2,  the 
Reynolds  number  is  defined  as  Rej  =  uRj / u,  where  Rj  is  the  exit  radius,  u  is  the  exit  flow 
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velocity  and  u  is  the  kinematic  viscosity  of  the  medium  used.  Whitehead  (1968)  used  the 
peak  pulsating  velocity,  while  the  mean  pulsating  velocity  was  selected  as  u  in  Sallet  and 
Widmayer’s  work  (1974)  and  the  present  study.  It  also  appears  that  there  is  no  significant 
effect  of  cross-flow  on  vortex  ring  formation  under  the  present  experimental  range. 

3.5  Dimensional  Analysis  for  Pulsating  Jet  Trajectories 


If  the  data  for  the  different  velocity  ratios  (U j/Uoo)  are  collapsed  onto  a  single  line, 
simple  power  law  relations  describing  the  pulsating  jet  trajectories  in  the  cross-flow  can  be 
easily  determined.  This  approach  is  composed  of  the  combination  of  dimensional  analysis 
and  fundamental  physical  arguments  about  the  characteristics  of  each  jet  flow.  First  of  all, 
an  appropriate  length  scale  associated  with  the  jet  flow  characteristics  needs  to  be  defined. 

In  the  present  study,  the  dependent  parameter  of  major  importance  is  the  penetration 
of  the  vortex  ring  in  the  vertical  direction.  It  can  be  described  by  the  following  function 

g{Uj,dj,Uoo,UT,i',Z,X,tJ,ld)^0  (7) 


where  /  is  a  pulsating  frequency.  Id  is  the  separation  distance  between  adjacent  vortex 
rings.  Dimensional  analysis  gives 


( VI  }U  Ujdj  Z  x\ 


(8) 


in  which  the  Strouhal  number  is  defined  as 


The  selection  of  Id  and  Ut  rather  than  dj  and  U j  for  the  definition  of  Strouhal  number 
is  that  the  vortex  ring  motion  is  more  dependent  on  translational  velocity  and  separation 
distance  than  jet  velocity  and  diameter. 

For  pulsed  jets,  the  Reynolds  number  based  on  the  pulsation  velocity  may  be  relevant 
in  the  formation  of  vortex  rings.  However,  since  no  efforts  have  been  made  on  the  Reynolds 
number  effect  in  this  study,  it  will  be  omitted  from  here  on. 
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Therefore,  the  penetration  can  be  described  by 


The  vertical  trajectory  of  a  vortex  ring  jet  is  expected  to  depend  on  the  Strouhal  number 
based  upon  translational  velocity  and  separation  distance  between  adjacent  vortex  rings 
and  the  momentum  flux  (U jfU^). 

In  order  to  correlate  the  pulsating  jet  flow  penetration  for  different  jet  free  stream 
velocity  ratios  and  various  frequencies,  a  modified  frequency-momentum  length  scale  is 
introduced  similar  to  the  steady  case 


where  k  is  a  similarity  parameter  to  be  determined  through  a  least  square  fit  of  the  data. 
The  data  of  Chassaing  et  al.  (1974)  was  fitted  using  a  regression  process.  In  this  calculation 
the  K  resulted  as  3.8.  The  result  in  Figure  9  indicates  a  good  correlation  of  all  the  various 
data.  This  correlation  is  being  considered  as  a  new  similarity  coordinate.  From  the  match 
between  experimental  results  and  similarity  analysis,  we  can  estimate  that  the  effects  of 
both  excitation  and  the  ratio  of  momentum  flux  (U j  on  the  trajectory  are  dominant 

for  the  pulsating  jet.  It  has  been  well  known  that  the  momentum  flux  ratio  (UjjU^)  is 
one  main  parameter  correlating  the  characteristics  of  the  steady  jet  in  a  cross-flow. 


3.6  Blooming  Rings 

An  interesting  observation  was  obtained  in  the  course  of  investigating  the  vortex  rings 
in  cross-flow.  At  very  low  cross-flow  speeds  {UjfUoo  >  3)  and  for  10  Hz  frequency  the 
vortex  ring  trajectory  was  not  stable  and  the  rings  were  found  to  frequently  move  in  various 
directions  into  the  cross-flow  and  into  the  upstream  direction,  in  a  lateral  direction  and  the 
various  orientations.  Figure  10.  The  resulting  flow  appeared  like  a  blooming  jet  obtained 
by  combined  axial  and  orbital  excitations  (Lee  and  Reynolds,  1984).  These  observations 
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confirmed  that  a  blooming  jet  can  be  generated  by  implementation  of  a  certain  geometry 
to  provide  passive  helical  excitation  coincident  with  a  resonance  mode  of  the  pipe  without 
active  orbital  excitation.  The  divergence  angle  is  a  reflection  of  the  amplitude  of  the  orbital 
(helical)  excitation.  The  divergence  angle  of  the  blooming  vortex  rings  was  smaller  than 
that  of  Lee  and  Reynolds  (1984).  The  mechanism  for  blooming  jets  and  bifurcating  jets  is 
discussed  by  Parekh  et  ah  (1988)  and  Lee  and  Reynolds  (1984).  The  prominent  mechanism 
for  these  particular  excited  jets  is  known  to  be  due  to  the  interaction  of  the  vortex  with 
the  initial  conditions. 

Anderson  (1955)  reported  in  a  paper  on  the  structure  and  velocity  of  the  periodic 
vortex  ring  pattern  of  a  Pfeifenton  (Pipe  Tone)  jet.  He  found  different  and  very  unusual 
vortex  flow  patterns,  which  showed  an  amazing  and  reproducible  splitting  of  the  jet  into 
two  distinct  circular  vortex  trains  (bifurcating  jet)  at  around  1150  cps.  In  his  paper 
he  described  the  phenomena  without  explanation  of  the  mechanism  —  this  supports  the 
notion  that  it  is  possible  to  obtain  bifurcating  and/or  blooming  jets  by  passive  excitation 
as  shown  in  this  experiment.  This  is  important  in  helping  to  understand  how  the  rings  are 
generated  in  such  excited  jets. 
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4.  Combustion  of  Pulsed  Fuel  Jets 


Based  on  the  experimental  and  numerical  results  obtained  (to  be  discussed  in  the  next 
section)  it  is  reasonably  clear  that,  if  the  flow  is  pulsed  such  that  vortical  rings  are  formed 
as  the  jet  is  issued  into  the  cross-flow,  then  increased  penetration  is  assured.  However,  two 
questions  still  persisted:  one.  Will  vortex  rings  also  be  formed  once  a  jet  is  pulsed  into  a 
supersonic  cross-flow?  and  two.  What  will  be  the  combustion  effect  on  the  rings?  To  find 
an  answer  to  these  questions  the  following  experiments  were  performed. 

Pulsed  jet  injection  into  a  supersonic  cross-flow  was  attempted  to  investigate  possi¬ 
ble  formation  of  vortex  rings.  An  extensive  search  for  a  reasonably  sized,  fast  response, 
commercially  available  valve  was  not  successful.  In  the  absence  of  an  off-the-shelf  valve 
to  pulse  the  flow,  a  small  free  piston  shock  tube  with  an  exit  diameter  of  1.25  cm  was 
designed  and  manufactured.  Figure  11.  The  tube  was  installed  in  the  tunnel  such  that 
its  main  axis  was  perpendicular  to  the  tunnel  stream  wise  direction.  High  speed  schlieren 
photography  with  a  spark  light  source  revealed  the  leading  shock  details,  wake  structures, 
and  the  jet  evolution  as  flow  emerged  from  the  open  end  of  the  shock  tube.  A  fast  speed 
camera  was  optically  coupled  with  the  shock  tube  to  record  the  injection  process.  The 
camera  took  eight  pictures  with  100  microseconds  delay  between  two  sequential  photos. 
The  starting  trigger  mechanism  for  the  camera  could  be  timed  for  different  initial  delays, 
this  was  necessary  because  of  the  desire  to  record  various  phases  of  the  jet  issued  into  the 
cross-flow. 

Figures  12,  13  and  14  show  the  sequence  of  a  pulsed  flow  out  of  the  tube  into  the 
tunnel  with  and  without  supersonic  flow.  A  delay  time  was  used  to  allow  the  jet  to  leave 
the  tube,  before  it  is  photographed.  In  these  figures  the  time  delay  for  the  camera  was 
about  500  to  650  microseconds,  each  frame  is  taken  in  20  microseconds,  and  each  frame  is 
separated  from  the  next  by  a  time  interval  of  100  microseconds.  As  the  flow  is  emerging 
from  the  jet  exit  (frame  2)  at  high  speeds,  a  shock  wave  is  formed.  However,  behind  the 
shock  wave  appears  to  be  a  turbulent  mushroom  shape  ring  vortex.  This  is  interpreted 
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to  be  a  supersonic  vortex  ring  which  we  have  intended  to  generate.  The  development  of 
this  single  pulsed  jet  flow  can  be  traced  in  sequential  photos  shown  in  this  figure  (frames 
2  through  8).  Figure  13  shows  a  similar  vortex  process  within  the  tunnel  with  cross-flow 
of  nominal  Mach  1.8  established.  In  Figure  13  the  initial  time  delay  is  slightly  more  and 
the  flow  is  leaving  the  jet  exit  in  the  first  frame.  As  expected,  an  oblique  shock  is  formed 
just  upstream  of  the  jet  exit  flow.  Even  though  the  pictures  are  not  so  clear,  a  supersonic 
vortex  ring  (mushroom)  is  visible  in  the  jet  flow.  The  ring  is  only  clearly  visible  in  the 
frames  1  and  2  of  the  Figure  13.  From  the  shock  structure  deformation,  the  approximate 
position  of  the  ring  can  be  approximately  identified.  Figure  14  shows  another  shot  similar 
as  the  one  in  Figure  13,  except  with  an  initial  delay  time  of  650  microseconds.  There 
appears  to  be  strong  similarities  between  the  behavior  of  vortex  rings  in  supersonic  flow 
and  that  of  low  speed  flows  considered  earlier. 

A  valve  capable  of  being  continuously  pulsed  at  very  high  rates,  orders  up  to  kilohertz, 
and  necessary  to  maintain  combustion  could  not  be  found  commercially.  An  effort  to 
design  an  in-house  valve  was  abandoned  after  it  was  realized  that  it  may  become  a  project 
of  uncontrolled  and  large  magnitude  by  itself.  There  could  be  no  comparison  made  with 
the  steady  jet  flow  for  the  penetration  into  the  cross-flow,  since  the  injection  mechanism 
would  not  lend  itself  to  be  used  as  a  steady  state  flow.  It  is  believed  that  sufficient  and 
accurate  data  for  steady  flow  injection  into  a  supersonic  cross-flow  exists  and  could  be  used 
for  comparison,  if  needed.  For  example,  these  would  be  consistent  with  the  previous  data 
on  Figure  9. 

Therefore,  an  effort  was  undertaken  to  study  the  combustion  of  a  pulsed  propane  fuel 
jet  in  ambient  air.  The  main  point  being  the  investigation  of  vortex  rings  formation  and 
motion  during  the  combustion  process.  A  setup  using  the  piezoelectric  valve,  (General 
Valve  model  9),  noted  earlier,  was  used  to  investigate  the  evolution  of  pulsed  propane 
(C3H8)  jets.  A  high  voltage  electric  arc,  between  two  thin  copper  needles.  Figure  15,  at 
the  exit  of  the  jet,  was  used  as  the  spark  to  initiate  combustion.  As  noted  earlier,  here  too. 
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the  intent  was  to  determine,  if  coherent  structures  such  as  vortex  rings  are  formed  during 
combustion  due  to  the  pulsation  of  the  fuel  jet  into  the  ambient  flow.  A  Kodak  fast  video 
recorder,  at  the  rate  of  one  thousand  frames  per  second  was  used  to  record  images  of  the 
pulsating  methane  gas  into  the  ambient  air.  As  observed  from  images  shown  in  Figure  16, 
one  can  only  cautiously  identify  the  burning  puffs  of  the  fuel  as  a  vortex  ring,  even  though 
periodically  there  are  clean  puffs  resembling  a  burning  turbulent  vortex  ring.  Results  of 
this  visual  step  could  be  tainted  because  the  mechanism  with  which  the  combustion  was 
initiated  was  two  dimensional  and  not  axisymmetric.  However,  there  is  clear  evidence  that 
under  certain  conditions,  vortex  rings  are  formed  and  that  the  motion  of  these  rings  in 
a  cross-flow  is  of  significant  interest.  This  is  relevant,  since  the  low  speed  study  revealed 
such  a  significant  increase  in  the  penetration  in  comparison  with  steady  jets. 
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5.  Numerical  Study 


The  main  purpose  of  this  numerical  study  is  to  better  understand  the  interaction 
between  vortex  rings  and  cross-flow,  and  to  explain  the  flow  visualization  experiments. 

Vortex  rings  in  cross-flow  represent  typical  three-dimensional,  unsteady  vortical  flow. 
Vortex  rings  are  usually  characterized  by  low  viscosity  and  concentrated  regions  of  vorticity 
embedded  in  an  otherwise  potential  flow.  This  fact  allows  one  to  discretize  the  compact 
regions  of  vorticity  into  an  assembly  of  vortices.  Vortex  methods  simulate  this  type  of 
flow  by  discretizing  and  tracking  the  discretization  in  a  Lagrangian  reference  frame.  A 
simplified  flow  simulation  was  carried  out  using  the  discrete  vortex  method.  That  is,  a 
vortex  ring  was  modeled  by  a  vortex  filament  with  a  number  of  nodes  in  the  shape  of 
a  torus.  A  discrete  number  of  nodes  which  lie  along  the  vortex  centerline  enable  one  to 
be  able  to  solve  the  deformation  problem.  These  points  are  located  along  a  given  initial 
configuration  of  the  ring,  along  with  the  initial  conditions  on  the  translational  velocity  of 
each  mode. 

Each  node  on  the  vortex  ring  moves  with  the  local  velocity.  Therefore,  once  the 
induced  velocities  are  known,  the  movement  of  each  node  during  a  small  time  step  can  be 
calculated  by  multiplying  the  local  induced  velocity  by  the  time  step.  For  this  calculation, 
a  second  order  Runge-Kutta  method  was  used  to  iterate  on  the  position  of  the  nodes  and 
therefore,  improved  the  accuracy  of  the  computations. 

The  original  code  was  written  and  developed  by  Parekh  et  al.  (1988)  based  upon  the 
discrete-vortex  method.  This  code  was  developed  to  investigate  the  effects  of  excitation 
frequency  and  amplitude  on  the  spreading  of  the  bifurcating  jet.  They  used  an  analytical 
function  to  describe  the  source  flow  and  discrete,  computational  vortex  elements  to  rep¬ 
resent  the  vortex  rings  formed  at  the  jet  exit.  This  function  was  modified  for  our  study 
to  investigate  the  effect  of  cross-flow  on  the  vortex  rings  depending  on  axial  frequency. 
We  computed  the  evolution  of  vortex  rings  in  the  cross-flow.  The  advantage  of  the  vortex 
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method  is  that  no  finite  mesh  or  grid  is  required,  and  this  Lagrangian  nature  makes  it 
unconditionally  stable  to  convection. 

The  viscosity  effect  was  not  initially  incorporated  in  this  numerical  experiment.  As  a 
result,  it  was  not  possible  to  obtain  the  characteristics  of  three-dimensional  ring  configu¬ 
ration,  including  shape  and  size.  In  order  to  complement  this  defect,  jet  source  flow  was 
introduced,  whose  flow  from  a  semi-infinite  cylindrical  vortex  sheet  is  similar  to  that  of  a 
viscous  jet  issuing  from  a  wall,  Chang  (1991). 


5.1  Formulation  and  Numerical  Scheme 


The  equations  which  govern  an  incompressible,  in  viscid  flow  are  Euler’s  equations; 

Du  du  _  ^  ^ 

_  =  =  (10) 

and  the  continuity  equation, 

V-ti  =  0  (11) 

expressing  the  conservation  of  momentum  and  mass,  respectively. 

This  computational  experiment  is  to  simulate  an  incompressible  flow  field  of  vortex 
rings  where  only  a  small  portion  is  rotational.  It  will  be  convenient  to  describe  the  above 
governing  equations  in  terms  of  vorticity  variable.  Equation  (10)  can  be  rewritten  in  terms 
of  the  vorticity  by  taking  the  curl. 


Du  du 


Dt 


-f  u  •  Vw  =  u  •  Vu 


dt 

a>  =  V  X  u 


(12) 

(13) 


Equation  (12)  is  called  the  vorticity  transport  equation  which  explains  that  the  vorticity 
moves  along  a  particle  path  while  it  is  being  tilted  and  stretched  with  the  evolving  strain 
field  (Vu).  Therefore,  the  term  [u  ■  Vu)  exists  only  in  a  three-dimensional  flow. 

The  vorticity  distribution  from  m  filaments  or  tubes  with  circulation  F  in  three- 
dimensional  flow  is  given  by 

dri{s) 

A  i  f  li 

1=1 


m  « 

=  X^Fi  /  7i[r,ri(s,t)] 


ds 


ds 


(14) 
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where  r{s,t)  describes  the  space  curve  of  the  filament  centerline  in  terms  of  a  parameter 
along  the  curve  (s).  The  distribution  function  ji  is  assumed  to  have  the  following  form 


(15) 


where  the  notation  of  Leonard  is  used  with  <7,-  (the  core  radius  of  filament  i). 


Once  the  vorticity  distribution  is  known,  the  velocity  induced  by  the  vorticity  in  a 
bounding  region  can  be  calculated  from  the  integration  of  Eqs.  (11)  and  (13).  We  can 
drive  using  vector  identity  and  V  •  u  =  0  that  (see  Batchelor^®) 


=  —V  X  LO 


(16) 


u  can  be  expressed  as 


(r  —  r')  X  U3{r')dr' 


r  —  r 


/|3 


+  V<^ 


(17) 


where  (j)  is  the  velocity  potential  related  to  the  homogeneous  solution  of  the  above  equation 
which  is  required  to  satisfy  the  boundary  conditions  as  follows: 


u  •  n 


=  0 


body  surface 


(18) 


lim  =u{r,t)  =  Uoo  (19) 

irHoo 


The  resultant  velocity  will  be  split  into  two  components  based  upon  the  uniqueness  of  the 
decomposition  of  a  vector  field. 


‘Ur  =  Usir,t)  +  ui 


(20) 


where  Us(r,t)  is  a  solenoidal  field  which  is  a  solution  of  Eq.  (16)  and  Uj  is  an  external 
irrotational  velocity.  The  induced  velocity  of  the  filament  at  point  r(s)  on  the  axis  is  given 
by  Biot-Savart  law 


Us{r,t) 


f{r-r')xu{r',t) 

|r  -r'P 


dV(r') 


(21) 
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With  the  thin  filament  approximation,  it  reduces  to 


/  r  f  [r{s,t)-r'{s\t)] 

Uair.t)  =  /  - 

Jc  \r{s,t) -r'{s', 


or 


(22) 


where  s  describes  the  coordinate  along  the  lines.  The  result  of  this  locally  induced  velocities 
which  in  the  limit  of  a  vortex  filament  becomes  infinite  yields  a  logarithmic  singularity  in 
the  self-induced  velocity  known  as  a  general  feature  of  a  curved  vortex  filament  in  the 
three-dimensional  flow.  This  causes  an  infinite  velocity  of  the  filament  and  also  causes 
instabilities.  In  order  to  eliminate  this  velocity  singularities,  numerous  cut  off  schemes 
have  been  introduced,  (see  Sarpkaya^^).  Another  difficulty  with  the  Biot-Savart  method 
will  be  the  number  of  operations  for  the  induced  velocity  calculation.  Because  CPU  time 
increases  considerably  as  the  number  of  vortex  elements  increase. 

The  modification  of  Biot-Savart  law  is  achieved  by  Rosenhead’s^^  method  which  was 
modified  by  Moore^^  later.  The  evolution  equation  for  a  three-dimensional,  space  curves 
=  l,2,...,iV)  for  N  vortex  filaments.  The  induced  velocity  at  a  point,  r,  is 
calculated  as: 

Us{r,t)  = 


[r{s,t)-rj{s',t)]  X 


,2 


+  (xafj 


(23) 


{\r{s,t)-rj\ 

where  acr^-  is  a  measure  of  the  vortical  flow  region,  (jj  is  the  core  radius  and  can  be  changed 
depending  on  the  distribution  of  vorticity.  A  vortex  ring  with  a  Gaussian  vorticity  distri¬ 
bution  can  be  simulated  by  a  single  vortex  filament  with  the  above  vorticity  distribution, 
and  then  the  self-induced  velocity  has  correct  value  for  a  =  0.413,  (Leonard^'*).  Therefore, 
this  value  was  used  for  our  numerical  simulations.  Here,  a  is  a  parameter  related  to  the 
fraction  of  circulation  within  the  radius  r  =  a  (Moore*^^). 

To  consider  the  influence  of  vortex  filament  i  on  filament  ji,  can  be  replaced  by 


oi{af  -f  crj)/2  as  follows: 
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(24) 


f  [rt(.s,t)  —  rj(6',t)]  x  ■’gy  ’  ^ 

j  [lri(s,i) -rjf +Q;'(cr? +crj)]  ^ 

As  mentioned  before,  a  =  0.413  satisfies  the  translational  velocity  of  a  vortex  ring  as  far 
as  ai  is  much  less  than  ring  radius. 

Since  a  vortex  line  always  consists  of  the  same  fluid  particle  and  moves  with  the  fluid 
in  an  incompressible  flow,  the  core  radius  a  is  kept  uniform  along  each  filament  but  varies 
with  time  to  conserve  a  constant  volume.  It  can  be  expressed  by 


(25) 


where  Ri{t)  is  the  instantaneous  length  of  the  vortex  ring  filament  i.  The  size  of  the  core 
radius  can  vary  with  the  local  stretching  of  the  ring.  The  total  circulation  assigned  to  a 
vortex  segment  at  each  simulation  maintains  constant. 

5.2  Calculation  Procedures 

The  dynamics  of  vortex  ring  in  cross-flow  is  characterized  by  three-dimensional  con¬ 
figuration  which  has  variations  with  location,  size  and  shape.  For  this  simulation,  a  set 
of  nodes  (vortex  elements)  are  incorporated  to  represent  each  vortex  ring.  The  spacing 
between  node  points  was  sufficiently  small  so  that  length  scales  of  interest  would  be  rep¬ 
resented  accurately. 

For  numerical  purposes,  it  is  assumed  that  the  core  radius  of  each  filament  is  much 
smaller  than  the  ring  radius.  Therefore,  the  distribution  of  vorticity  is  almost  Gaussian  in 
magnitude  and  the  vortex  lines  are  parallel  at  each  cross  section  normal  to  the  filament 
tangent  vector. 

Since  each  vortex  element  moves  with  the  local  fluid  velocity,  the  self-induced  velocity 
field  obtained  from  Eq.  (24)  can  be  integrated  in  time  to  find  the  subsequent  location  of  the 
vortex  elements.  Again  using  Eq.  (24)  the  velocity  field  at  the  next  time  can  be  computed 
and  this  procedure  repeated. 
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The  total  velocity  induced  by  a  multi-segment  filament  is  obtained  by  summing  the 
induced  velocity  of  each  filament.  The  detail  calculations  for  the  induced  velocity  by  all 
the  vortex  filaments  are  given  in  Parekh.^^ 


The  first  vortex  ring  is  generated  at  the  origin  with  the  initialization  of  frequency 
and  time  integration  parameters.  Then,  subsequent  vortex  rings  are  generated  at  given 
distances  depending  on  axial  frequency.  The  circulation  of  each  filament  has  the  same 
magnitude  and  is  determined  by  circulation  conservation  laws.  The  fiux  of  vorticity  in  the 
boundary  layer  (sheet)  is  given  by  17^/2.  It  is  also  assumed  that  the  vorticity  convected 
from  the  cylindrical  sheet  equals  the  vorticity  convected  by  the  discrete  filaments.  The 
conservation  constraints  show 


r 

At 


(26) 


in  which  P  is  the  circulation  of  each  ring  filament,  7  is  the  circulation  per  unit  length  of 
the  cylindrical  vortex  sheet,  and  At  is  the  time  interval  between  generation  of  vortex  ring 
filaments. 

At  each  time  step,  the  total  velocity  at  each  node  due  to  the  combined  effects  of  the 
vortex  segments,  jet  source  flow  and  external  flow  is  calculated.  Both  the  spatial  and  the 
temporal  integration  are  carried  out  with  the  second-order,  Runge— Kutta  method,  with 
the  node  points  being  advanced  in  time 

I'=x(()+<i(()|  (27) 

x{t  -|-  At)  =  x(t)  +  u'(t)6t  (28) 


in  which  x{t)  and  u{t)  indicate  a  node  position  and  velocity  at  time  t,  and  6t  is  the  time 
increment  between  time  steps  which  is  typically  an  order  of  magnitude  smaller  than  At. 


5.3  Simulation  of  Jet  Source  Flow 

The  formulation  presented  below  follows  the  work  of  Parekh.'*®  In  order  to  simulate 
the  jet  source  flow,  the  Biot-Savart  integral  is  applied  to  a  semi-infinite,  cylindrical  vortex 


-27- 


sheet  of  finite  thickness.  The  jet  velocity  function  can  be  derived  to  define  the  mean  jet 
flow.  The  simple  form  of  equation  is  given  by 

Ur(or,?)  =  (29) 


where 


Uz(or,2)  = 


7  TT  1(2  —  ovh) 


27ror/i  2 


+  ovh 


T  +  t; 


+  {H  -  ox)T{n,vi) 


1/  =  1  +  or^  +  +  acr^, 


mi  =  1  — 


(1  -^) 


(1+^)’ 

ovh  =  1  +  or^  +  , 


m  ^  r 

Jo 

T(u,rj)=  r 
Jo 


Jo  (1  —  u  cos  6)^/1  —  fj  cos  $ 

The  detail  derivation  of  Eqs.  (29)  and  (30)  and  integration  scheme  based  upon  the 
fourth-order  Newton-Cotes  method  are  described  in  Ref.  45. 

The  semi-infinite  sheet  of  vorticity  extends  from  — oo  to  origin  as  shown  in  Fig.  19. 
Its  axis  represents  the  jet  centerline,  and  the  end  of  the  sheet  represents  the  jet  exit. 
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Here,  it  is  assumed  that  the  cylindrical  sheet  effects  the  motion  of  the  filaments  but  there 
is  no  effect  of  the  filaments  on  the  cylindrical  sheet.  The  trajectory  of  each  filament  is 
determined  by  the  total  velocities  induced  by  each  filament,  jet  source  function  and  by  a 
uniform  cross-flow.  The  sheet  maintains  cylindrical  configuration  and  it  is  not  moved  due 
to  the  filaments  since  the  cylindrical  sheet  corresponds  to  a  mean  flow  whose  centerline  is 
determined  only  by  the  position  of  a  physical  nozzle. 

5.4  Numerical  Results 

Numerical  results  are  presented  to  investigate  the  interaction  between  vortex  rings  and 
cross-flows.  The  remarkable  features  of  the  trajectories  and  motion  of  vortex  rings  with 
different  effects  of  cross-flow,  circulation  strength,  and  Strouhal  number  are  illustrated  as 
follows. 

In  order  to  simulate  the  dynamics  of  vortex  rings  in  a  cross-flow,  discrete  vortex 
filaments  are  introduced  to  represent  the  vortex  rings.  Each  vortex  filament  consists  of  32 
nodes  or  vortex  elements,  and  each  vortex  segment  moves  under  the  action  of  the  induced 
velocity  of  the  other  segments  and  the  cross-flow.  Since  the  variation  of  core  radius  has 
a  very  small  effect  on  the  motion  of  vortex  rings  in  a  cross-flow,  the  uniform  core  radius 
(a  =  0.1)  was  selected  in  this  calculation.  The  non-dimensional  time  step  for  the  evolution 
of  the  vortex  ring  was  typically  below  8t  =  0.08.  Twenty-one  different  cases  were  tested  to 
simulate  as  closely  as  possible  the  experimental  results.  The  stream  wise  cross-flow  is  in 
the  X-direction,  Z  is  the  vertical  direction  measured  from  the  jet  center.  An  initial  vortex 
ring  is  generated  at  the  origin,  and  then  the  next  ring,  begins  to  evolve  under  the  effect  of 
the  induced  velocity. 

In  Figure  17,  the  numerical  results  show  some  typical  and  similar  tilted  vortex  rings 
having  the  same  tilted  direction  as  in  the  experiment.  The  downstream  portion  of  the 
vortex  ring  moves  upward  under  the  effect  of  non  uniform  interactions  between  the  vortex 
ring  and  the  cross-flow.  At  each  time,  the  side  and  end  views  are  represented  in  Figure  18. 
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The  large  value  of  time  was  used  to  observe  the  deformation  of  large  numbers  of  vortex 
rings  on  each  simulation.  In  the  following  subsections,  the  effects  of  pulsations  through 
the  change  of  Strouhal  number  (frequency),  circulation  strength  and  cross-flow  including 
constantly  uniform,  uniform  shear  and  non-uniform  flows  will  be  described  in  terms  of 
ring  trajectories  and  motions.  The  sequences  of  events  during  vortex  ring  evolution  are 
shown  in  Figures  19  and  20  as  a  function  of  frequency.  The  data  in  each  figure  was  given 
a  different  magnitude  of  uniform  velocity.  As  the  frequency  increases,  the  spatial  distance 
between  neighboring  vortex  rings  becomes  smaller.  Also,  the  penetration  at  relatively  high 
frequency  is  reduced  due  to  the  strong  interactions  between  adjacent  vortex  rings.  This 
observation  was  revealed  in  the  experimental  work  described  above.  Therefore,  it  is  seen 
that  the  spatial  distance  between  vortex  rings  is  determined  by  the  axial  pulsation  (through 
Strouhal  number).  Naturally,  the  total  number  of  vortex  rings  generated  are  increased 
within  the  fixed  coordinate  frame  as  the  axial  frequency  increases.  Also,  the  increase 
of  uniform  velocity  makes  the  vortex  rings  more  deflected  and  consequently  reduces  the 
penetration  with  the  same  pulsation  conditions.  This  causes  the  total  number  of  rings  to  be 
decreased  within  the  fixed  frame.  The  averaged  tilting  angle  in  this  numerical  calculation 
is  in  good  agreement  with  the  measured  value  (about  3°  to  4°)  except  for  the  case  of 
St  =  0.35  and  Uoo  =  0.125. 


The  initial  circulation  strength  of  the  vortex  ring  has  a  significant  effect  on  the  tra¬ 
jectories  of  rings  in  the  cross- flow.  As  illustrated  in  Figures  21  cind  22,  the  penetration 
of  vortex  rings  into  the  cross-flow  are  considerably  deeper  with  the  increase  of  circulation 
strength.  The  circulation  effect  was  investigated  for  two  different  ranges  of  Strouhal  num¬ 
ber.  At  relatively  high  Strouhal  number  {St  =  0.35),  as  the  circulation  increases,  a  tangle 
of  vortex  rings  is  caused  due  to  the  strong  interaction  of  two  neighboring  vortex  rings  as 
presented  in  Figure  23.  In  Figure  24  {St  =  0.25),  as  expected,  a  considerable  trend  of 
increase  of  penetration  was  observed  without  entanglement  by  the  interaction  of  adjacent 
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vortex  rings.  Also,  the  tilted  configuration  of  the  ring  was  not  considerably  changed  as 
long  as  the  relatively  low  pulsation  is  maintained. 

We  also  investigated  the  effect  of  cross-flow  on  the  vortex  ring,  by  just  adding  dif¬ 
ferent  kinds  of  velocity  profiles.  It  was  discussed  by  Aref  and  Flinchem  (1984)  that  the 
incorporation  of  the  background  flow  is  an  approximation  which  neglects  the  effect  of  the 
vortex  ring  on  the  background  flow.  The  case  of  uniform  shear  flow  was  considered  first 
with  U{Z)  =  mZ  where  m  is  0.125.  Next,  we  impose  a  non  uniform  shear  flow  with  the 
exponential  velocity  profile  as  follows, 

UiZ)  =  n{l-exp{-Z))  (12) 

With  n  =  0.125,  this  profile  has  the  same  shear  rate  at  Z  =  20  as  in  the  uniform  shear 
flow  case.  This  indicates  that  it  has  larger  shear  rates  than  the  uniform  shear  case  below 
Z  —  20.  The  vortex  ring  evolution  is  shown  in  Figures  25  and  26.  The  uniform  shear 
flow  has  more  penetration  due  to  the  low  shear  rates,  however,  at  relatively  high  Strouhal 
numbers  it  tends  to  easily  entangle  with  strong  interactions  between  adjacent  vortex  rings. 
The  exponential  velocity  profile  shows  more  or  less  close  trajectories  and  motion  of  vortex 
rings  compared  to  the  constantly  uniform  cross-flow  case. 

Two  important  results  are  illustrated  from  this  numerical  simulation.  One  is  the 
various  penetrations  of  vortex  rings  due  to  the  effect  of  pulsations,  circulations,  and  cross- 
flows.  The  other  is  the  tilted  configuration  of  vortex  rings  in  the  cross-flow.  This  numerical 
simulation  confirmed  that  the  tilting  mechanism  of  the  vortex  ring  is  due  to  the  combined 
effect  of  cross-flow  and  the  viscosity  of  the  jet  source  flow.  As  shown  in  the  experimental 
results,  the  basic  circular  shape  is  maintained  within  the  fixed  coordinate  frame.  In  order 
to  make  more  meaningful  comparisons  between  numerical  and  experimental  results,  it  is 
necessary  to  match  basic  parameters  and  initial  conditions  for  the  generation  of  vortex 
rings  accurately.  That  is,  at  the  same  time  t,  the  position  and  induced  velocity  of  the 
vortex  ring  should  be  calculated  to  reasonably  predict  the  subsequent  behavior  of  the  ring. 
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6,  Conclusions 


The  following  conclusions  have  been  drawn  from  this  research,  based  on  the  exper¬ 
imental  observations  including  measurements  and  computational  results.  As  a  result  of 
full  pulsation  of  the  flow  of  a  steady  jet,  vortex  rings  were  formed.  At  low  pulsing  fre¬ 
quencies  the  vortex  rings  penetrated  into  the  cross-flow  much  deeper  than  for  high  pulsing 
frequencies  or  for  the  steady  jet.  Interactions  between  vortex  rings  became  dominant  at 
low  jet  to  cross-flow  velocity  ratios  and  for  high  frequencies.  No  appreciable  interaction 
was  observed  between  neighboring  vortex  rings  below  pulsing  frequencies  of  4  Hz  for  a  jet 
to  cross-flow  velocity  ratio  of  three.  The  vortex  rings  at  low  frequencies  moved  as  unit 
blobs  of  fluid  with  little  change  in  their  general  shape.  Analyses  showed  that  pulsation 
frequency  and  momentum  flux  effects  were  two  factors  dominating  dynamics  of  pulsating 
jet  flow.  Penetration  trajectories  of  the  pulsed  jets  (vortex  rings)  were  characterized  by 
a  power  law  for  each  velocity  ratio  and  pulsing  frequency.  Pulsation  of  the  steady  jet, 
therefore,  resulted  in  increased  penetration  and  incresed  mixing. 

The  velocity  field  measured  within  vortex  rings,  generated  by  pulsing  at  1  Hz  for  a 
jet  to  cross-flow  velocity  ratio  of  3.0,  revealed  fundamental  characteristic  features  of  single 
vortex  rings  in  cross-flow.  Measurements  provided  information  that  vortex  rings  are  fully- 
formed  at  Z  =  3dj  which  is  in  good  agreement  with  the  observations  by  others,  in  still  air, 
as  well  as  computational  results.  Circulation  measured  by  the  hot-film  probe  was  close  to 
that  predicted  by  a  simple  ‘slug  model’.  At  Z  =  Zdj,  9.2%  of  the  ideal  circulation  was  lost 
from  the  ring  while  24.2%  was  lost  at  lOd^. 

Computational  results  supported  and  agreed  with  the  increased  penetration  of  the 
pulsed  jet  and  the  tilting  of  individual  rings,  as  observed  in  the  experiments.  As  the  fre¬ 
quency  (Strouhal  number)  increased,  the  spatial  distance  between  the  vortex  rings  became 
smaller,  and  the  penetration  was  reduced  due  to  strong  interactions  between  adjacent 
vortex  rings  as  observed  from  the  experiments. 
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Table  1 

Quantitative  data  from  hot-film  measurement 


z 

Zdj 

7d,- 

lOdj 

Ring  radius  i2,cm 

1.080 

1.250 

1.397 

Translational  velocity 

9.635 

9.162 

8.471 

Radius  ratio  a/R 

0.294 

0.279 

0.268 

Measured  circulation  rni,cm^/^cc 

38.526 

35.183 

33.631 

Core  circulation  Fccm^/sec 

25.896 

23.326 

21.615 

Ideal  circulation  ri,CTn^/sec 

39.879 

39.879 

39.879 

Mean  pulsating  velocity  u,  cm/sec 

12.630 

12.630 

12.630 

Reynolds  number  Re  =  ^ 

3970 

3970 

3970 

Tilting  angle 

2.6* 

3.2® 

3.3® 

ln[}^)  -  ^  in  inviscid  equation 
Non-dimensional  velocity  U 

3.054 

3.106 

3.146 

3.610 

4.411 

4.919 

X-coordinate  of  ring  center 

0.40<fj  ■ 

1.24<ij 

l.ndj 

Table  2 

Comparative  data  on  vortex  ring  formation 


Ring  formation  distance 

Re„ 

Medium 

Whitehead  (1968) 

2.0d, 

533 

Computation 

Sallet  et  al.  (1974) 

2,Uj  . 

8010  -  12250 

Experiment  in  air 

This  study 

3.0d, 

800 

Experiment  in  water 
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Figure  1.  Experimental  setup  of  water  tunnel  for  continuous  pulsation. 


(a)  1  Hz  (b)  10  Hz 


Figure  2.  Typical  pulsation  profiles  at  the  exit. 


Figure  3.  Schematic  diagram  for  thermal  anemometer  and  data  acquisition  system 


(a)  Pulsed  at  1  Hz 


(b)  Pulsed  at  2  Hz 


(c)  Pulsed  at  10  Hz 


(d)  Steady  jet 


(e)  Top  view  of  a  ring  at  1  Hz 


Figure  4. 


41 


Figure  4.  Typical  photographs  of  vortex  rings  in  cross-flow  for  various  pulsation  and  flow 
conditions  for  Uj/Uoo  =  1.5,  Red-  =  1870 

(a)  Pulsed  at  1  Hz,  (b)  Pulsed  at  2  Hz,  (c)  Pulsed  at  10  Hz,  (d)  Steady  Jet, 

(e)  Top  view  of  a  ring,  jet  is  pulsed  at  1  Hz. 

(f)  A  single  ring  moving  into  the  cross-flow,  note  the  upstream  motion  caused  by 
the  tilting  process 
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(s)  -0.696  it 


(k)  Xm^o,mR 


(i)  (0 


Figure  7.  Typical  velocity  flow-field  within  a  ring  at  .Z  =  3dj 

(a)  through  (i)  are  schematically  explained  in  Figure  4.7(j) 
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Figure  8.  Sequence  of  traces  showing  fully-formed  vortex  ring 


Figure  10.  Blooming  jet  observed  without  any  active  control  of  the  exit  geometry. 
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Tunnel  Flow,  M>  1 


Figure  11.  Schematic  of  the  shock  tube  used  to  create  a  pulsed  flow  into  a  high  speed 
cross-flow 
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Figure  12.  Sequence  of  a  pulsed  jet  comming  out  of  the  shock  tube  into  still  air.  Each 
frame  is  taken  in  20  microseconds  with  100  microseconds  between  sequential  frames. 
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Figure  13.  Sequence  of  a  pulsed  jet  coming  out  of  the  shock  tube  into  a  1.80  Mach  cross 
flow  of  air.  Initial  delay  time  of  550  microsecond.  Each  frame  is  taken  in  20  microseconds 
with  100  microseconds  between  sequential  frames. 
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Figure  14.  Sequence  of  a  pulsed  jet  coming  out  of  the  shock  tube  into  a  1.80  Mach  cross 
flow  of  air.  Initial  delay  time  of  650  microsecond.  Each  frame  is  taken  in  20  microseconds 
with  100  microseconds  between  sequential  frames. 
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60mm 


Figure  16.  Combusting  images  of  a  pulsed  fuel  jet  into  ambient  air  at  20  Hz.  Burning 
vortex  rings  are  apparent. 
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Figure  22.  Effect  of  circulation  on  motion  and  trajectory  at  St  =  0.25, 


- 

Cross-flow 

cross- flow 

direction 

velocity  profile 

_ 1 _ 1 _ 1 _ mlnm - 

-0.0 

a  X 

.0 

(a)  Uniform  flow 

-5.0  20.0  ^ 


(b)  Mean  shear  flow 


(c)  Exponential  flow 


Figure  23.  Effect  of  cross-flow  on  motion  and  trajectory  at  St  =  0.35, 
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Figure  24.  Effect  of  cross-flow  on  motion  and  trajectory  at  St  =  0.25 
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Figure  25.  Time  evolution  of  jet  pulses,  issued  in  a  uniform  cross-flow, 
and  the  formation  of  vortex  rings. 


(b)  t  =  20.0 


Figure  26.  Time  evolution  of  jet  pulses,  issued  in  a  shear  cross-flow,  and  the  formation  of 
vortex  rings. 


